------------------------------------------------------------------------
-- ddSubtract.decTest -- decDouble subtraction                        --
-- Copyright (c) IBM Corporation, 1981, 2008.  All rights reserved.   --
------------------------------------------------------------------------
-- Please see the document "General Decimal Arithmetic Testcases"     --
-- at http://www2.hursley.ibm.com/decimal for the description of      --
-- these testcases.                                                   --
--                                                                    --
-- These testcases are experimental ('beta' versions), and they       --
-- may contain errors.  They are offered on an as-is basis.  In       --
-- particular, achieving the same results as the tests here is not    --
-- a guarantee that an implementation complies with any Standard      --
-- or specification.  The tests are not exhaustive.                   --
--                                                                    --
-- Please send comments, suggestions, and corrections to the author:  --
--   Mike Cowlishaw, IBM Fellow                                       --
--   IBM UK, PO Box 31, Birmingham Road, Warwick CV34 5JL, UK         --
--   mfc@uk.ibm.com                                                   --
------------------------------------------------------------------------
version: 2.59

-- This set of tests are for decDoubles only; all arguments are
-- representable in a decDouble
precision:   16
maxExponent: 384
minExponent: -383
extended:    1
clamp:       1
rounding:    half_even

-- [first group are 'quick confidence check']
ddsub001 subtract  0   0  -> '0'
ddsub002 subtract  1   1  -> '0'
ddsub003 subtract  1   2  -> '-1'
ddsub004 subtract  2   1  -> '1'
ddsub005 subtract  2   2  -> '0'
ddsub006 subtract  3   2  -> '1'
ddsub007 subtract  2   3  -> '-1'

ddsub011 subtract -0   0  -> '-0'
ddsub012 subtract -1   1  -> '-2'
ddsub013 subtract -1   2  -> '-3'
ddsub014 subtract -2   1  -> '-3'
ddsub015 subtract -2   2  -> '-4'
ddsub016 subtract -3   2  -> '-5'
ddsub017 subtract -2   3  -> '-5'

ddsub021 subtract  0  -0  -> '0'
ddsub022 subtract  1  -1  -> '2'
ddsub023 subtract  1  -2  -> '3'
ddsub024 subtract  2  -1  -> '3'
ddsub025 subtract  2  -2  -> '4'
ddsub026 subtract  3  -2  -> '5'
ddsub027 subtract  2  -3  -> '5'

ddsub030 subtract  11  1  -> 10
ddsub031 subtract  10  1  ->  9
ddsub032 subtract  9   1  ->  8
ddsub033 subtract  1   1  ->  0
ddsub034 subtract  0   1  -> -1
ddsub035 subtract -1   1  -> -2
ddsub036 subtract -9   1  -> -10
ddsub037 subtract -10  1  -> -11
ddsub038 subtract -11  1  -> -12

ddsub040 subtract '5.75' '3.3'  -> '2.45'
ddsub041 subtract '5'    '-3'   -> '8'
ddsub042 subtract '-5'   '-3'   -> '-2'
ddsub043 subtract '-7'   '2.5'  -> '-9.5'
ddsub044 subtract '0.7'  '0.3'  -> '0.4'
ddsub045 subtract '1.3'  '0.3'  -> '1.0'
ddsub046 subtract '1.25' '1.25' -> '0.00'

ddsub050 subtract '1.23456789'    '1.00000000' -> '0.23456789'
ddsub051 subtract '1.23456789'    '1.00000089' -> '0.23456700'

ddsub060 subtract '70'    '10000e+16' -> '-1.000000000000000E+20' Inexact Rounded
ddsub061 subtract '700'    '10000e+16' -> '-1.000000000000000E+20' Inexact Rounded
ddsub062 subtract '7000'    '10000e+16' -> '-9.999999999999999E+19' Inexact Rounded
ddsub063 subtract '70000'    '10000e+16' -> '-9.999999999999993E+19' Rounded
ddsub064 subtract '700000'    '10000e+16' -> '-9.999999999999930E+19' Rounded
  -- symmetry:
ddsub065 subtract '10000e+16'    '70' -> '1.000000000000000E+20' Inexact Rounded
ddsub066 subtract '10000e+16'    '700' -> '1.000000000000000E+20' Inexact Rounded
ddsub067 subtract '10000e+16'    '7000' -> '9.999999999999999E+19' Inexact Rounded
ddsub068 subtract '10000e+16'    '70000' -> '9.999999999999993E+19' Rounded
ddsub069 subtract '10000e+16'    '700000' -> '9.999999999999930E+19' Rounded

  -- some of the next group are really constructor tests
ddsub090 subtract '00.0'    '0.0'  -> '0.0'
ddsub091 subtract '00.0'    '0.00' -> '0.00'
ddsub092 subtract '0.00'    '00.0' -> '0.00'
ddsub093 subtract '00.0'    '0.00' -> '0.00'
ddsub094 subtract '0.00'    '00.0' -> '0.00'
ddsub095 subtract '3'    '.3'   -> '2.7'
ddsub096 subtract '3.'   '.3'   -> '2.7'
ddsub097 subtract '3.0'  '.3'   -> '2.7'
ddsub098 subtract '3.00' '.3'   -> '2.70'
ddsub099 subtract '3'    '3'    -> '0'
ddsub100 subtract '3'    '+3'   -> '0'
ddsub101 subtract '3'    '-3'   -> '6'
ddsub102 subtract '3'    '0.3'  -> '2.7'
ddsub103 subtract '3.'   '0.3'  -> '2.7'
ddsub104 subtract '3.0'  '0.3'  -> '2.7'
ddsub105 subtract '3.00' '0.3'  -> '2.70'
ddsub106 subtract '3'    '3.0'  -> '0.0'
ddsub107 subtract '3'    '+3.0' -> '0.0'
ddsub108 subtract '3'    '-3.0' -> '6.0'

-- the above all from add; massaged and extended.  Now some new ones...
-- [particularly important for comparisons]
-- NB: -xE-8 below were non-exponents pre-ANSI X3-274, and -1E-7 or 0E-7
-- with input rounding.
ddsub120 subtract  '10.23456784'    '10.23456789'  -> '-5E-8'
ddsub121 subtract  '10.23456785'    '10.23456789'  -> '-4E-8'
ddsub122 subtract  '10.23456786'    '10.23456789'  -> '-3E-8'
ddsub123 subtract  '10.23456787'    '10.23456789'  -> '-2E-8'
ddsub124 subtract  '10.23456788'    '10.23456789'  -> '-1E-8'
ddsub125 subtract  '10.23456789'    '10.23456789'  -> '0E-8'
ddsub126 subtract  '10.23456790'    '10.23456789'  -> '1E-8'
ddsub127 subtract  '10.23456791'    '10.23456789'  -> '2E-8'
ddsub128 subtract  '10.23456792'    '10.23456789'  -> '3E-8'
ddsub129 subtract  '10.23456793'    '10.23456789'  -> '4E-8'
ddsub130 subtract  '10.23456794'    '10.23456789'  -> '5E-8'
ddsub131 subtract  '10.23456781'    '10.23456786'  -> '-5E-8'
ddsub132 subtract  '10.23456782'    '10.23456786'  -> '-4E-8'
ddsub133 subtract  '10.23456783'    '10.23456786'  -> '-3E-8'
ddsub134 subtract  '10.23456784'    '10.23456786'  -> '-2E-8'
ddsub135 subtract  '10.23456785'    '10.23456786'  -> '-1E-8'
ddsub136 subtract  '10.23456786'    '10.23456786'  -> '0E-8'
ddsub137 subtract  '10.23456787'    '10.23456786'  -> '1E-8'
ddsub138 subtract  '10.23456788'    '10.23456786'  -> '2E-8'
ddsub139 subtract  '10.23456789'    '10.23456786'  -> '3E-8'
ddsub140 subtract  '10.23456790'    '10.23456786'  -> '4E-8'
ddsub141 subtract  '10.23456791'    '10.23456786'  -> '5E-8'
ddsub142 subtract  '1'              '0.999999999'  -> '1E-9'
ddsub143 subtract  '0.999999999'    '1'            -> '-1E-9'
ddsub144 subtract  '-10.23456780'   '-10.23456786' -> '6E-8'
ddsub145 subtract  '-10.23456790'   '-10.23456786' -> '-4E-8'
ddsub146 subtract  '-10.23456791'   '-10.23456786' -> '-5E-8'

-- additional scaled arithmetic tests [0.97 problem]
ddsub160 subtract '0'     '.1'      -> '-0.1'
ddsub161 subtract '00'    '.97983'  -> '-0.97983'
ddsub162 subtract '0'     '.9'      -> '-0.9'
ddsub163 subtract '0'     '0.102'   -> '-0.102'
ddsub164 subtract '0'     '.4'      -> '-0.4'
ddsub165 subtract '0'     '.307'    -> '-0.307'
ddsub166 subtract '0'     '.43822'  -> '-0.43822'
ddsub167 subtract '0'     '.911'    -> '-0.911'
ddsub168 subtract '.0'    '.02'     -> '-0.02'
ddsub169 subtract '00'    '.392'    -> '-0.392'
ddsub170 subtract '0'     '.26'     -> '-0.26'
ddsub171 subtract '0'     '0.51'    -> '-0.51'
ddsub172 subtract '0'     '.2234'   -> '-0.2234'
ddsub173 subtract '0'     '.2'      -> '-0.2'
ddsub174 subtract '.0'    '.0008'   -> '-0.0008'
-- 0. on left
ddsub180 subtract '0.0'     '-.1'      -> '0.1'
ddsub181 subtract '0.00'    '-.97983'  -> '0.97983'
ddsub182 subtract '0.0'     '-.9'      -> '0.9'
ddsub183 subtract '0.0'     '-0.102'   -> '0.102'
ddsub184 subtract '0.0'     '-.4'      -> '0.4'
ddsub185 subtract '0.0'     '-.307'    -> '0.307'
ddsub186 subtract '0.0'     '-.43822'  -> '0.43822'
ddsub187 subtract '0.0'     '-.911'    -> '0.911'
ddsub188 subtract '0.0'     '-.02'     -> '0.02'
ddsub189 subtract '0.00'    '-.392'    -> '0.392'
ddsub190 subtract '0.0'     '-.26'     -> '0.26'
ddsub191 subtract '0.0'     '-0.51'    -> '0.51'
ddsub192 subtract '0.0'     '-.2234'   -> '0.2234'
ddsub193 subtract '0.0'     '-.2'      -> '0.2'
ddsub194 subtract '0.0'     '-.0008'   -> '0.0008'
-- negatives of same
ddsub200 subtract '0'     '-.1'      -> '0.1'
ddsub201 subtract '00'    '-.97983'  -> '0.97983'
ddsub202 subtract '0'     '-.9'      -> '0.9'
ddsub203 subtract '0'     '-0.102'   -> '0.102'
ddsub204 subtract '0'     '-.4'      -> '0.4'
ddsub205 subtract '0'     '-.307'    -> '0.307'
ddsub206 subtract '0'     '-.43822'  -> '0.43822'
ddsub207 subtract '0'     '-.911'    -> '0.911'
ddsub208 subtract '.0'    '-.02'     -> '0.02'
ddsub209 subtract '00'    '-.392'    -> '0.392'
ddsub210 subtract '0'     '-.26'     -> '0.26'
ddsub211 subtract '0'     '-0.51'    -> '0.51'
ddsub212 subtract '0'     '-.2234'   -> '0.2234'
ddsub213 subtract '0'     '-.2'      -> '0.2'
ddsub214 subtract '.0'    '-.0008'   -> '0.0008'

-- more fixed, LHS swaps [really the same as testcases under add]
ddsub220 subtract '-56267E-12' 0  -> '-5.6267E-8'
ddsub221 subtract '-56267E-11' 0  -> '-5.6267E-7'
ddsub222 subtract '-56267E-10' 0  -> '-0.0000056267'
ddsub223 subtract '-56267E-9'  0  -> '-0.000056267'
ddsub224 subtract '-56267E-8'  0  -> '-0.00056267'
ddsub225 subtract '-56267E-7'  0  -> '-0.0056267'
ddsub226 subtract '-56267E-6'  0  -> '-0.056267'
ddsub227 subtract '-56267E-5'  0  -> '-0.56267'
ddsub228 subtract '-56267E-2'  0  -> '-562.67'
ddsub229 subtract '-56267E-1'  0  -> '-5626.7'
ddsub230 subtract '-56267E-0'  0  -> '-56267'
-- symmetry ...
ddsub240 subtract 0 '-56267E-12'  -> '5.6267E-8'
ddsub241 subtract 0 '-56267E-11'  -> '5.6267E-7'
ddsub242 subtract 0 '-56267E-10'  -> '0.0000056267'
ddsub243 subtract 0 '-56267E-9'   -> '0.000056267'
ddsub244 subtract 0 '-56267E-8'   -> '0.00056267'
ddsub245 subtract 0 '-56267E-7'   -> '0.0056267'
ddsub246 subtract 0 '-56267E-6'   -> '0.056267'
ddsub247 subtract 0 '-56267E-5'   -> '0.56267'
ddsub248 subtract 0 '-56267E-2'   -> '562.67'
ddsub249 subtract 0 '-56267E-1'   -> '5626.7'
ddsub250 subtract 0 '-56267E-0'   -> '56267'

-- now some more from the 'new' add
ddsub301 subtract '1.23456789'  '1.00000000' -> '0.23456789'
ddsub302 subtract '1.23456789'  '1.00000011' -> '0.23456778'

-- some carrying effects
ddsub321 subtract '0.9998'  '0.0000' -> '0.9998'
ddsub322 subtract '0.9998'  '0.0001' -> '0.9997'
ddsub323 subtract '0.9998'  '0.0002' -> '0.9996'
ddsub324 subtract '0.9998'  '0.0003' -> '0.9995'
ddsub325 subtract '0.9998'  '-0.0000' -> '0.9998'
ddsub326 subtract '0.9998'  '-0.0001' -> '0.9999'
ddsub327 subtract '0.9998'  '-0.0002' -> '1.0000'
ddsub328 subtract '0.9998'  '-0.0003' -> '1.0001'

-- internal boundaries
ddsub346 subtract '10000e+9'  '7'   -> '9999999999993'
ddsub347 subtract '10000e+9'  '70'   -> '9999999999930'
ddsub348 subtract '10000e+9'  '700'   -> '9999999999300'
ddsub349 subtract '10000e+9'  '7000'   -> '9999999993000'
ddsub350 subtract '10000e+9'  '70000'   -> '9999999930000'
ddsub351 subtract '10000e+9'  '700000'   -> '9999999300000'
ddsub352 subtract '7' '10000e+9'   -> '-9999999999993'
ddsub353 subtract '70' '10000e+9'   -> '-9999999999930'
ddsub354 subtract '700' '10000e+9'   -> '-9999999999300'
ddsub355 subtract '7000' '10000e+9'   -> '-9999999993000'
ddsub356 subtract '70000' '10000e+9'   -> '-9999999930000'
ddsub357 subtract '700000' '10000e+9'   -> '-9999999300000'

-- zero preservation
ddsub361 subtract 1 '0.0001' -> '0.9999'
ddsub362 subtract 1 '0.00001' -> '0.99999'
ddsub363 subtract 1 '0.000001' -> '0.999999'
ddsub364 subtract 1 '0.0000000000000001' -> '0.9999999999999999'
ddsub365 subtract 1 '0.00000000000000001' -> '1.000000000000000' Inexact Rounded
ddsub366 subtract 1 '0.000000000000000001' -> '1.000000000000000' Inexact Rounded

-- some funny zeros [in case of bad signum]
ddsub370 subtract 1  0  -> 1
ddsub371 subtract 1 0.  -> 1
ddsub372 subtract 1  .0 -> 1.0
ddsub373 subtract 1 0.0 -> 1.0
ddsub374 subtract  0  1 -> -1
ddsub375 subtract 0.  1 -> -1
ddsub376 subtract  .0 1 -> -1.0
ddsub377 subtract 0.0 1 -> -1.0

-- leading 0 digit before round
ddsub910 subtract -103519362 -51897955.3 -> -51621406.7
ddsub911 subtract 159579.444 89827.5229 -> 69751.9211

ddsub920 subtract 333.0000000123456 33.00000001234566 -> 299.9999999999999 Inexact Rounded
ddsub921 subtract 333.0000000123456 33.00000001234565 -> 300.0000000000000 Inexact Rounded
ddsub922 subtract 133.0000000123456 33.00000001234565 ->  99.99999999999995
ddsub923 subtract 133.0000000123456 33.00000001234564 ->  99.99999999999996
ddsub924 subtract 133.0000000123456 33.00000001234540 -> 100.0000000000002 Rounded
ddsub925 subtract 133.0000000123456 43.00000001234560 ->  90.00000000000000
ddsub926 subtract 133.0000000123456 43.00000001234561 ->  89.99999999999999
ddsub927 subtract 133.0000000123456 43.00000001234566 ->  89.99999999999994
ddsub928 subtract 101.0000000123456 91.00000001234566 ->   9.99999999999994
ddsub929 subtract 101.0000000123456 99.00000001234566 ->   1.99999999999994

-- more LHS swaps [were fixed]
ddsub390 subtract '-56267E-10'   0 ->  '-0.0000056267'
ddsub391 subtract '-56267E-6'    0 ->  '-0.056267'
ddsub392 subtract '-56267E-5'    0 ->  '-0.56267'
ddsub393 subtract '-56267E-4'    0 ->  '-5.6267'
ddsub394 subtract '-56267E-3'    0 ->  '-56.267'
ddsub395 subtract '-56267E-2'    0 ->  '-562.67'
ddsub396 subtract '-56267E-1'    0 ->  '-5626.7'
ddsub397 subtract '-56267E-0'    0 ->  '-56267'
ddsub398 subtract '-5E-10'       0 ->  '-5E-10'
ddsub399 subtract '-5E-7'        0 ->  '-5E-7'
ddsub400 subtract '-5E-6'        0 ->  '-0.000005'
ddsub401 subtract '-5E-5'        0 ->  '-0.00005'
ddsub402 subtract '-5E-4'        0 ->  '-0.0005'
ddsub403 subtract '-5E-1'        0 ->  '-0.5'
ddsub404 subtract '-5E0'         0 ->  '-5'
ddsub405 subtract '-5E1'         0 ->  '-50'
ddsub406 subtract '-5E5'         0 ->  '-500000'
ddsub407 subtract '-5E15'        0 ->  '-5000000000000000'
ddsub408 subtract '-5E16'        0 ->  '-5.000000000000000E+16'  Rounded
ddsub409 subtract '-5E17'        0 ->  '-5.000000000000000E+17'  Rounded
ddsub410 subtract '-5E18'        0 ->  '-5.000000000000000E+18'  Rounded
ddsub411 subtract '-5E100'       0 ->  '-5.000000000000000E+100' Rounded

-- more RHS swaps [were fixed]
ddsub420 subtract 0  '-56267E-10' ->  '0.0000056267'
ddsub421 subtract 0  '-56267E-6'  ->  '0.056267'
ddsub422 subtract 0  '-56267E-5'  ->  '0.56267'
ddsub423 subtract 0  '-56267E-4'  ->  '5.6267'
ddsub424 subtract 0  '-56267E-3'  ->  '56.267'
ddsub425 subtract 0  '-56267E-2'  ->  '562.67'
ddsub426 subtract 0  '-56267E-1'  ->  '5626.7'
ddsub427 subtract 0  '-56267E-0'  ->  '56267'
ddsub428 subtract 0  '-5E-10'     ->  '5E-10'
ddsub429 subtract 0  '-5E-7'      ->  '5E-7'
ddsub430 subtract 0  '-5E-6'      ->  '0.000005'
ddsub431 subtract 0  '-5E-5'      ->  '0.00005'
ddsub432 subtract 0  '-5E-4'      ->  '0.0005'
ddsub433 subtract 0  '-5E-1'      ->  '0.5'
ddsub434 subtract 0  '-5E0'       ->  '5'
ddsub435 subtract 0  '-5E1'       ->  '50'
ddsub436 subtract 0  '-5E5'       ->  '500000'
ddsub437 subtract 0  '-5E15'      ->  '5000000000000000'
ddsub438 subtract 0  '-5E16'      ->  '5.000000000000000E+16'   Rounded
ddsub439 subtract 0  '-5E17'      ->  '5.000000000000000E+17'   Rounded
ddsub440 subtract 0  '-5E18'      ->  '5.000000000000000E+18'   Rounded
ddsub441 subtract 0  '-5E100'     ->  '5.000000000000000E+100'  Rounded


-- try borderline precision, with carries, etc.
ddsub461 subtract '1E+16' '1'        -> '9999999999999999'
ddsub462 subtract '1E+12' '-1.111'   -> '1000000000001.111'
ddsub463 subtract '1.111'  '-1E+12'  -> '1000000000001.111'
ddsub464 subtract '-1'    '-1E+16'   -> '9999999999999999'
ddsub465 subtract '7E+15' '1'        -> '6999999999999999'
ddsub466 subtract '7E+12' '-1.111'   -> '7000000000001.111'
ddsub467 subtract '1.111'  '-7E+12'  -> '7000000000001.111'
ddsub468 subtract '-1'    '-7E+15'   -> '6999999999999999'

--                  1234567890123456       1234567890123456      1 23456789012345
ddsub470 subtract '0.4444444444444444'  '-0.5555555555555563' -> '1.000000000000001' Inexact Rounded
ddsub471 subtract '0.4444444444444444'  '-0.5555555555555562' -> '1.000000000000001' Inexact Rounded
ddsub472 subtract '0.4444444444444444'  '-0.5555555555555561' -> '1.000000000000000' Inexact Rounded
ddsub473 subtract '0.4444444444444444'  '-0.5555555555555560' -> '1.000000000000000' Inexact Rounded
ddsub474 subtract '0.4444444444444444'  '-0.5555555555555559' -> '1.000000000000000' Inexact Rounded
ddsub475 subtract '0.4444444444444444'  '-0.5555555555555558' -> '1.000000000000000' Inexact Rounded
ddsub476 subtract '0.4444444444444444'  '-0.5555555555555557' -> '1.000000000000000' Inexact Rounded
ddsub477 subtract '0.4444444444444444'  '-0.5555555555555556' -> '1.000000000000000' Rounded
ddsub478 subtract '0.4444444444444444'  '-0.5555555555555555' -> '0.9999999999999999'
ddsub479 subtract '0.4444444444444444'  '-0.5555555555555554' -> '0.9999999999999998'
ddsub480 subtract '0.4444444444444444'  '-0.5555555555555553' -> '0.9999999999999997'
ddsub481 subtract '0.4444444444444444'  '-0.5555555555555552' -> '0.9999999999999996'
ddsub482 subtract '0.4444444444444444'  '-0.5555555555555551' -> '0.9999999999999995'
ddsub483 subtract '0.4444444444444444'  '-0.5555555555555550' -> '0.9999999999999994'

-- and some more, including residue effects and different roundings
rounding: half_up
ddsub500 subtract '1231234567456789' 0             -> '1231234567456789'
ddsub501 subtract '1231234567456789' 0.000000001   -> '1231234567456789' Inexact Rounded
ddsub502 subtract '1231234567456789' 0.000001      -> '1231234567456789' Inexact Rounded
ddsub503 subtract '1231234567456789' 0.1           -> '1231234567456789' Inexact Rounded
ddsub504 subtract '1231234567456789' 0.4           -> '1231234567456789' Inexact Rounded
ddsub505 subtract '1231234567456789' 0.49          -> '1231234567456789' Inexact Rounded
ddsub506 subtract '1231234567456789' 0.499999      -> '1231234567456789' Inexact Rounded
ddsub507 subtract '1231234567456789' 0.499999999   -> '1231234567456789' Inexact Rounded
ddsub508 subtract '1231234567456789' 0.5           -> '1231234567456789' Inexact Rounded
ddsub509 subtract '1231234567456789' 0.500000001   -> '1231234567456788' Inexact Rounded
ddsub510 subtract '1231234567456789' 0.500001      -> '1231234567456788' Inexact Rounded
ddsub511 subtract '1231234567456789' 0.51          -> '1231234567456788' Inexact Rounded
ddsub512 subtract '1231234567456789' 0.6           -> '1231234567456788' Inexact Rounded
ddsub513 subtract '1231234567456789' 0.9           -> '1231234567456788' Inexact Rounded
ddsub514 subtract '1231234567456789' 0.99999       -> '1231234567456788' Inexact Rounded
ddsub515 subtract '1231234567456789' 0.999999999   -> '1231234567456788' Inexact Rounded
ddsub516 subtract '1231234567456789' 1             -> '1231234567456788'
ddsub517 subtract '1231234567456789' 1.000000001   -> '1231234567456788' Inexact Rounded
ddsub518 subtract '1231234567456789' 1.00001       -> '1231234567456788' Inexact Rounded
ddsub519 subtract '1231234567456789' 1.1           -> '1231234567456788' Inexact Rounded

rounding: half_even
ddsub520 subtract '1231234567456789' 0             -> '1231234567456789'
ddsub521 subtract '1231234567456789' 0.000000001   -> '1231234567456789' Inexact Rounded
ddsub522 subtract '1231234567456789' 0.000001      -> '1231234567456789' Inexact Rounded
ddsub523 subtract '1231234567456789' 0.1           -> '1231234567456789' Inexact Rounded
ddsub524 subtract '1231234567456789' 0.4           -> '1231234567456789' Inexact Rounded
ddsub525 subtract '1231234567456789' 0.49          -> '1231234567456789' Inexact Rounded
ddsub526 subtract '1231234567456789' 0.499999      -> '1231234567456789' Inexact Rounded
ddsub527 subtract '1231234567456789' 0.499999999   -> '1231234567456789' Inexact Rounded
ddsub528 subtract '1231234567456789' 0.5           -> '1231234567456788' Inexact Rounded
ddsub529 subtract '1231234567456789' 0.500000001   -> '1231234567456788' Inexact Rounded
ddsub530 subtract '1231234567456789' 0.500001      -> '1231234567456788' Inexact Rounded
ddsub531 subtract '1231234567456789' 0.51          -> '1231234567456788' Inexact Rounded
ddsub532 subtract '1231234567456789' 0.6           -> '1231234567456788' Inexact Rounded
ddsub533 subtract '1231234567456789' 0.9           -> '1231234567456788' Inexact Rounded
ddsub534 subtract '1231234567456789' 0.99999       -> '1231234567456788' Inexact Rounded
ddsub535 subtract '1231234567456789' 0.999999999   -> '1231234567456788' Inexact Rounded
ddsub536 subtract '1231234567456789' 1             -> '1231234567456788'
ddsub537 subtract '1231234567456789' 1.00000001    -> '1231234567456788' Inexact Rounded
ddsub538 subtract '1231234567456789' 1.00001       -> '1231234567456788' Inexact Rounded
ddsub539 subtract '1231234567456789' 1.1           -> '1231234567456788' Inexact Rounded
-- critical few with even bottom digit...
ddsub540 subtract '1231234567456788' 0.499999999   -> '1231234567456788' Inexact Rounded
ddsub541 subtract '1231234567456788' 0.5           -> '1231234567456788' Inexact Rounded
ddsub542 subtract '1231234567456788' 0.500000001   -> '1231234567456787' Inexact Rounded

rounding: down
ddsub550 subtract '1231234567456789' 0             -> '1231234567456789'
ddsub551 subtract '1231234567456789' 0.000000001   -> '1231234567456788' Inexact Rounded
ddsub552 subtract '1231234567456789' 0.000001      -> '1231234567456788' Inexact Rounded
ddsub553 subtract '1231234567456789' 0.1           -> '1231234567456788' Inexact Rounded
ddsub554 subtract '1231234567456789' 0.4           -> '1231234567456788' Inexact Rounded
ddsub555 subtract '1231234567456789' 0.49          -> '1231234567456788' Inexact Rounded
ddsub556 subtract '1231234567456789' 0.499999      -> '1231234567456788' Inexact Rounded
ddsub557 subtract '1231234567456789' 0.499999999   -> '1231234567456788' Inexact Rounded
ddsub558 subtract '1231234567456789' 0.5           -> '1231234567456788' Inexact Rounded
ddsub559 subtract '1231234567456789' 0.500000001   -> '1231234567456788' Inexact Rounded
ddsub560 subtract '1231234567456789' 0.500001      -> '1231234567456788' Inexact Rounded
ddsub561 subtract '1231234567456789' 0.51          -> '1231234567456788' Inexact Rounded
ddsub562 subtract '1231234567456789' 0.6           -> '1231234567456788' Inexact Rounded
ddsub563 subtract '1231234567456789' 0.9           -> '1231234567456788' Inexact Rounded
ddsub564 subtract '1231234567456789' 0.99999       -> '1231234567456788' Inexact Rounded
ddsub565 subtract '1231234567456789' 0.999999999   -> '1231234567456788' Inexact Rounded
ddsub566 subtract '1231234567456789' 1             -> '1231234567456788'
ddsub567 subtract '1231234567456789' 1.00000001    -> '1231234567456787' Inexact Rounded
ddsub568 subtract '1231234567456789' 1.00001       -> '1231234567456787' Inexact Rounded
ddsub569 subtract '1231234567456789' 1.1           -> '1231234567456787' Inexact Rounded

-- symmetry...
rounding: half_up
ddsub600 subtract 0             '1231234567456789' -> '-1231234567456789'
ddsub601 subtract 0.000000001   '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub602 subtract 0.000001      '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub603 subtract 0.1           '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub604 subtract 0.4           '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub605 subtract 0.49          '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub606 subtract 0.499999      '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub607 subtract 0.499999999   '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub608 subtract 0.5           '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub609 subtract 0.500000001   '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub610 subtract 0.500001      '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub611 subtract 0.51          '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub612 subtract 0.6           '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub613 subtract 0.9           '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub614 subtract 0.99999       '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub615 subtract 0.999999999   '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub616 subtract 1             '1231234567456789' -> '-1231234567456788'
ddsub617 subtract 1.000000001   '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub618 subtract 1.00001       '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub619 subtract 1.1           '1231234567456789' -> '-1231234567456788' Inexact Rounded

rounding: half_even
ddsub620 subtract 0             '1231234567456789' -> '-1231234567456789'
ddsub621 subtract 0.000000001   '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub622 subtract 0.000001      '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub623 subtract 0.1           '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub624 subtract 0.4           '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub625 subtract 0.49          '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub626 subtract 0.499999      '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub627 subtract 0.499999999   '1231234567456789' -> '-1231234567456789' Inexact Rounded
ddsub628 subtract 0.5           '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub629 subtract 0.500000001   '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub630 subtract 0.500001      '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub631 subtract 0.51          '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub632 subtract 0.6           '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub633 subtract 0.9           '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub634 subtract 0.99999       '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub635 subtract 0.999999999   '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub636 subtract 1             '1231234567456789' -> '-1231234567456788'
ddsub637 subtract 1.00000001    '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub638 subtract 1.00001       '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub639 subtract 1.1           '1231234567456789' -> '-1231234567456788' Inexact Rounded
-- critical few with even bottom digit...
ddsub640 subtract 0.499999999   '1231234567456788' -> '-1231234567456788' Inexact Rounded
ddsub641 subtract 0.5           '1231234567456788' -> '-1231234567456788' Inexact Rounded
ddsub642 subtract 0.500000001   '1231234567456788' -> '-1231234567456787' Inexact Rounded

rounding: down
ddsub650 subtract 0             '1231234567456789' -> '-1231234567456789'
ddsub651 subtract 0.000000001   '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub652 subtract 0.000001      '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub653 subtract 0.1           '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub654 subtract 0.4           '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub655 subtract 0.49          '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub656 subtract 0.499999      '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub657 subtract 0.499999999   '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub658 subtract 0.5           '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub659 subtract 0.500000001   '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub660 subtract 0.500001      '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub661 subtract 0.51          '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub662 subtract 0.6           '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub663 subtract 0.9           '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub664 subtract 0.99999       '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub665 subtract 0.999999999   '1231234567456789' -> '-1231234567456788' Inexact Rounded
ddsub666 subtract 1             '1231234567456789' -> '-1231234567456788'
ddsub667 subtract 1.00000001    '1231234567456789' -> '-1231234567456787' Inexact Rounded
ddsub668 subtract 1.00001       '1231234567456789' -> '-1231234567456787' Inexact Rounded
ddsub669 subtract 1.1           '1231234567456789' -> '-1231234567456787' Inexact Rounded


-- lots of leading zeros in intermediate result, and showing effects of
-- input rounding would have affected the following
rounding: half_up
ddsub670 subtract '1234567456789' '1234567456788.1' -> 0.9
ddsub671 subtract '1234567456789' '1234567456788.9' -> 0.1
ddsub672 subtract '1234567456789' '1234567456789.1' -> -0.1
ddsub673 subtract '1234567456789' '1234567456789.5' -> -0.5
ddsub674 subtract '1234567456789' '1234567456789.9' -> -0.9

rounding: half_even
ddsub680 subtract '1234567456789' '1234567456788.1' -> 0.9
ddsub681 subtract '1234567456789' '1234567456788.9' -> 0.1
ddsub682 subtract '1234567456789' '1234567456789.1' -> -0.1
ddsub683 subtract '1234567456789' '1234567456789.5' -> -0.5
ddsub684 subtract '1234567456789' '1234567456789.9' -> -0.9

ddsub685 subtract '1234567456788' '1234567456787.1' -> 0.9
ddsub686 subtract '1234567456788' '1234567456787.9' -> 0.1
ddsub687 subtract '1234567456788' '1234567456788.1' -> -0.1
ddsub688 subtract '1234567456788' '1234567456788.5' -> -0.5
ddsub689 subtract '1234567456788' '1234567456788.9' -> -0.9

rounding: down
ddsub690 subtract '1234567456789' '1234567456788.1' -> 0.9
ddsub691 subtract '1234567456789' '1234567456788.9' -> 0.1
ddsub692 subtract '1234567456789' '1234567456789.1' -> -0.1
ddsub693 subtract '1234567456789' '1234567456789.5' -> -0.5
ddsub694 subtract '1234567456789' '1234567456789.9' -> -0.9

-- Specials
ddsub780 subtract -Inf   Inf   -> -Infinity
ddsub781 subtract -Inf   1000  -> -Infinity
ddsub782 subtract -Inf   1     -> -Infinity
ddsub783 subtract -Inf  -0     -> -Infinity
ddsub784 subtract -Inf  -1     -> -Infinity
ddsub785 subtract -Inf  -1000  -> -Infinity
ddsub787 subtract -1000  Inf   -> -Infinity
ddsub788 subtract -Inf   Inf   -> -Infinity
ddsub789 subtract -1     Inf   -> -Infinity
ddsub790 subtract  0     Inf   -> -Infinity
ddsub791 subtract  1     Inf   -> -Infinity
ddsub792 subtract  1000  Inf   -> -Infinity

ddsub800 subtract  Inf   Inf   ->  NaN  Invalid_operation
ddsub801 subtract  Inf   1000  ->  Infinity
ddsub802 subtract  Inf   1     ->  Infinity
ddsub803 subtract  Inf   0     ->  Infinity
ddsub804 subtract  Inf  -0     ->  Infinity
ddsub805 subtract  Inf  -1     ->  Infinity
ddsub806 subtract  Inf  -1000  ->  Infinity
ddsub807 subtract  Inf  -Inf   ->  Infinity
ddsub808 subtract -1000 -Inf   ->  Infinity
ddsub809 subtract -Inf  -Inf   ->  NaN  Invalid_operation
ddsub810 subtract -1    -Inf   ->  Infinity
ddsub811 subtract -0    -Inf   ->  Infinity
ddsub812 subtract  0    -Inf   ->  Infinity
ddsub813 subtract  1    -Inf   ->  Infinity
ddsub814 subtract  1000 -Inf   ->  Infinity
ddsub815 subtract  Inf  -Inf   ->  Infinity

ddsub821 subtract  NaN   Inf   ->  NaN
ddsub822 subtract -NaN   1000  -> -NaN
ddsub823 subtract  NaN   1     ->  NaN
ddsub824 subtract  NaN   0     ->  NaN
ddsub825 subtract  NaN  -0     ->  NaN
ddsub826 subtract  NaN  -1     ->  NaN
ddsub827 subtract  NaN  -1000  ->  NaN
ddsub828 subtract  NaN  -Inf   ->  NaN
ddsub829 subtract -NaN   NaN   -> -NaN
ddsub830 subtract -Inf   NaN   ->  NaN
ddsub831 subtract -1000  NaN   ->  NaN
ddsub832 subtract -1     NaN   ->  NaN
ddsub833 subtract -0     NaN   ->  NaN
ddsub834 subtract  0     NaN   ->  NaN
ddsub835 subtract  1     NaN   ->  NaN
ddsub836 subtract  1000 -NaN   -> -NaN
ddsub837 subtract  Inf   NaN   ->  NaN

ddsub841 subtract  sNaN  Inf   ->  NaN  Invalid_operation
ddsub842 subtract -sNaN  1000  -> -NaN  Invalid_operation
ddsub843 subtract  sNaN  1     ->  NaN  Invalid_operation
ddsub844 subtract  sNaN  0     ->  NaN  Invalid_operation
ddsub845 subtract  sNaN -0     ->  NaN  Invalid_operation
ddsub846 subtract  sNaN -1     ->  NaN  Invalid_operation
ddsub847 subtract  sNaN -1000  ->  NaN  Invalid_operation
ddsub848 subtract  sNaN  NaN   ->  NaN  Invalid_operation
ddsub849 subtract  sNaN sNaN   ->  NaN  Invalid_operation
ddsub850 subtract  NaN  sNaN   ->  NaN  Invalid_operation
ddsub851 subtract -Inf -sNaN   -> -NaN  Invalid_operation
ddsub852 subtract -1000 sNaN   ->  NaN  Invalid_operation
ddsub853 subtract -1    sNaN   ->  NaN  Invalid_operation
ddsub854 subtract -0    sNaN   ->  NaN  Invalid_operation
ddsub855 subtract  0    sNaN   ->  NaN  Invalid_operation
ddsub856 subtract  1    sNaN   ->  NaN  Invalid_operation
ddsub857 subtract  1000 sNaN   ->  NaN  Invalid_operation
ddsub858 subtract  Inf  sNaN   ->  NaN  Invalid_operation
ddsub859 subtract  NaN  sNaN   ->  NaN  Invalid_operation

-- propagating NaNs
ddsub861 subtract  NaN01   -Inf     ->  NaN1
ddsub862 subtract -NaN02   -1000    -> -NaN2
ddsub863 subtract  NaN03    1000    ->  NaN3
ddsub864 subtract  NaN04    Inf     ->  NaN4
ddsub865 subtract  NaN05    NaN61   ->  NaN5
ddsub866 subtract -Inf     -NaN71   -> -NaN71
ddsub867 subtract -1000     NaN81   ->  NaN81
ddsub868 subtract  1000     NaN91   ->  NaN91
ddsub869 subtract  Inf      NaN101  ->  NaN101
ddsub871 subtract  sNaN011  -Inf    ->  NaN11  Invalid_operation
ddsub872 subtract  sNaN012  -1000   ->  NaN12  Invalid_operation
ddsub873 subtract -sNaN013   1000   -> -NaN13  Invalid_operation
ddsub874 subtract  sNaN014   NaN171 ->  NaN14  Invalid_operation
ddsub875 subtract  sNaN015  sNaN181 ->  NaN15  Invalid_operation
ddsub876 subtract  NaN016   sNaN191 ->  NaN191 Invalid_operation
ddsub877 subtract -Inf      sNaN201 ->  NaN201 Invalid_operation
ddsub878 subtract -1000     sNaN211 ->  NaN211 Invalid_operation
ddsub879 subtract  1000    -sNaN221 -> -NaN221 Invalid_operation
ddsub880 subtract  Inf      sNaN231 ->  NaN231 Invalid_operation
ddsub881 subtract  NaN025   sNaN241 ->  NaN241 Invalid_operation

-- edge case spills
ddsub901 subtract  2.E-3  1.002  -> -1.000
ddsub902 subtract  2.0E-3  1.002  -> -1.0000
ddsub903 subtract  2.00E-3  1.0020  -> -1.00000
ddsub904 subtract  2.000E-3  1.00200  -> -1.000000
ddsub905 subtract  2.0000E-3  1.002000  -> -1.0000000
ddsub906 subtract  2.00000E-3  1.0020000  -> -1.00000000
ddsub907 subtract  2.000000E-3  1.00200000  -> -1.000000000
ddsub908 subtract  2.0000000E-3  1.002000000  -> -1.0000000000

-- subnormals and overflows covered under Add

-- Null tests
ddsub9990 subtract 10  # -> NaN Invalid_operation
ddsub9991 subtract  # 10 -> NaN Invalid_operation
